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A GaAs/AlGaAs based two-qubit quantum device that allows the controlled generation and 
straightforward detection of entanglement by measuring a stationary current-voltage characteris- 
tic is proposed. We have developed a two-particle Green's function method of open systems and 
calculate the properties of three-dimensional interacting entangled systems non-perturbatively. We 
present concrete device designs and detailed, charge self-consistent predictions. One of the qubits 
is an all-electric Mach-Zehnder interferometer that consists of two electrostatically defined quan- 
tum wires with coupling windows, whereas the second qubit is an electrostatically defined double 
quantum dot located in a second two-dimensional electron gas beneath the quantum wires. We find 
that the entanglement of the device can be controlled externally by tuning the tunneling coupling 
between the two quantum dots. 
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I. INTRODUCTION 

Semiconductor based concepts for quantum informa- 
tion processing promise a high degree of scalability. How- 
ever, the implementation of even a single qubit proves to 
be more difficult in semiconductors than in alternative 
approaches. I'^i^iii^i^i The reasons lie in the short deco- 
herence times and strong interactions between elemen- 
tary excitations in solids4 While a 7-qubit NMR quan- 
tum computer was realized already in 2001,^ the first 
two-qubit quantum operation in semiconductors has been 
demonstrated with spin qubits only recently i^°i^^i^^ 

At liquid helium temperatures and below, electrons in 
a GaAs/AlGaAs two-dimensional electron gas (2DEG) 
show mean free paths as well as phase-relaxation lengths 
of the order of 10 — 20 /im which are remarkably long 
distancesJ^ Indeed, Bertoni et ali^ and lonicioiu et 
al.^^ have proposed a scheme for quantum computa- 
tion in semiconductors that exploits these long coher- 
ence lengths for electrons propagating through quantum 
wires (QWRs). In this scheme, a single electron prop- 
agates through two parallel QWRs, that represent the 
qubit states |0) and respectively. A single-qubit ro- 
tation gate can be realized by an electronic directional 
coupler— that may consist of a small window in the 
barrier between the two QWRs and is able to transfer 
the wave packet from one channel to the other. Single- 
qubit structures have been studied theoretically by sev- 
eral authorsi^iiSii^!^ who showed how the dimensions of 
the coupling window can be utilized to tailor the transfer 
process. This approach allows one, at least in principle, 
to apply multiple quantum gate operations to a qubit 
without exceeding the relevant coherence lengths. 

Universal quantum computation requires not only one- 
qubit but also two-qubit gates® that are responsible for 
the creation of entanglement between the qubits. To this 
end, two pairs of QWRs may be brought close to one 



another so that the electrons in the QWRs can interact 
with one another in a controlled way via their Coulomb 
repulsion. This concept has been analyzed both for wave 
packet922i2^iHi2^ as well as for stationary state a^°i^®'^'' . 
The latter investigations used simple models that provide 
important proofs-of-principle, but do not provide quan- 
titative predictions of realizable device structures that 
exhibit quantum gate operations. 

While two-qubit devices have only been addressed 
theoretically in the literature, one-qubit systems have 
also stimulated a number of experiments in both 
latera l^^i^^ as well as vertical^SiSi nanostructures. How- 
ever, quantum entanglement can clearly not be exposed 
in single-qubit systems as entanglement expresses the 
non-separability of a multi-qubit state. Therefore, the 
following step towards the realization of a prototype 
quantum gate based on coupled QWRs consists of a two- 
qubit device. Also, we should keep in mind that conceiv- 
able two-qubit experiments need realistic estimations of 
an entanglement witness that is experimentally accessible 
and a clear signature of the correct coherent functioning 
of the directional coupler. 

The aim of the present work is to theoretically study 
realistic single- and two-qubit quantum gates, by propos- 
ing conceptually simple and experimentally realizable 
novel semiconductor devices for quantum information 
processing. The system proposed is based on ballistic 
GaAs/AlGaAs QWRs and allows for the controlled gen- 
eration and detection of entanglement between an all- 
electric Mach-Zehnder interferometer and an electrostat- 
ically defined single-electron double quantum dot. The 
Mach-Zehnder interferometer is realized by two electro- 
statically defined QWRs that are connected by two cou- 
pling windows. We model the electrons that are propa- 
gating through the interferometer by stationary scatter- 
ing states, and our qubit state will be defined accordingly. 
In fact, the original proposal of Bertoni et at involves 
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the injection of electron wave packets into the pair of 
QWRs. However, we chose here a time-independent ap- 
proach since, (i) it is equivalent to the time-dependent 
one when the spatial dimension of the wave packet tends 
to be larger than the device (indeed, this is the case in 
our devices due to the small source/drain biases leading 
to a well defined kinetic energy of the electron), (ii) the 
structures that implement the quantum gates operate in 
the same way, (iii) it allows us to compute realistic es- 
timations of the charge-self-consistent ballistic I-V char- 
acteristics of the device. 

In our calculations, we include the detailed charge- 
self consistent three-dimensional device geometry, mate- 
rial composition, doping profile and bias voltage. Im- 
portantly, we demonstrate that the all-electric Mach- 
Zehnder interferometer can function as a fully control- 
lable single-qubit gate for experimentally attainable pa- 
rameters. In addition, we have developed a Green's func- 
tion method for the quantitative analysis of the entangled 
Mach-Zehnder and double quantum dot system that in- 
cludes the Coulomb interaction between the two qubits 
non-perturbatively. In order to gain better qualitative 
insight into the numerical results, we interpret them in 
terms of an analytical model that reproduces the com- 
putational results qualitatively. Taking both results to- 
gether, we are able to show that the degree of entangle- 
ment can be related to the DC I-V characteristics of the 
interferometer and that the Mach-Zehnder double quan- 
tum dot device can be employed as an externally control- 
lable two-qubit gate. 

This paper is organized as follows. In Sec. |TT1 we in- 
troduce a novel method for the quantum mechanical cal- 
culation of the ballistic current through an interacting 
two-particle system. This includes the calculation of the 
electronic structure and the determination of the bal- 
listic transmission characteristics of the realistic three- 
dimensional nanostructure. In Sec. IIIIl we discuss nu- 
merical details of the method. In Sec. IIVI we introduce 
a simplified analytical scheme able to gather the essen- 
tial features of the two-qubit device. The results obtained 
with this approach will be compared, in the following sec- 
tion, with the calculations of the realistic device. Sec. |V] 
focuses on the results and the discussion of the ballis- 
tic current through the proposed GaAs/AlGaAs single- 
and two-qubit QWR devices. In the same section, the 
degree of entanglement between the two qubits is eval- 
uated. The paper concludes with final remarks and a 
summary in Sec. I VII 



II. METHOD 

In this section, we present a novel method for the cal- 
culation of the stationary ballistic current of two inter- 
acting particles in realistic three-dimensional nanostruc- 
tures. We restrict ourselves to the case where only one of 
the two particles (say, particle 1) contributes to the cur- 
rent while the second particle (particle 2) is assumed to 



be bound. Importantly, this method accounts for many- 
particle effects non-perturbatively which is an important 
prerequisite of studying entanglement and, consequently, 
predicting solid-state based quantum information pro- 
cessing. The presently developed scheme generalizes the 
contact block reduction (CBR) method for single-particle 
ballistic currents that we have developed previously] '^^i'^^ 
In analogy to the CBR method, we proceed in two steps 
that we will describe in this section. First, we deter- 
mine the equilibrium electronic structure of a closed two- 
particle system, obtained by substituting the leads of our 
device with closed boundary conditions. Then, in the 
second step, we calculate the current carrying states of 
the open two-particle system. In fact, the closed calcula- 
tion only provides a convenient basis for the computation 
of the scattering states. For the sake of efficiency and 
without loss of generality, we perform the former calcu- 
lation in the Hartree approximation, by solving two cou- 
pled Schrodinger and Poisson equations. Then, we use 
single-particle product states as the basis for the scatter- 
ing states of the fully interacting two-particle device. 



A. Electronic structure of particle 1 

We consider a single-particle, single-band, effective 
mass Schrodinger equation for this electron, and include 
the electrostatic Hartree potential in the absence of the 
second particle. The closed system is represented by 
a Hermitian Hamiltonian with von Neumann boundary 
conditions for the wave functions at the contacts (van- 
ishing normal derivative) i^i^ The electronic structure 
of the closed device is calculated self-consistently. To 
this end, we iteratively solve the Schrodinger equation 
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and the non-linear Poisson equation 



Ve,. (x) £oV0 (x) 



(1) 



(2) 



until the electrostatic potential (f) and the total charge 
density p have reached convergence. The effective mass 
TO*, the conduction band offset E^, and the relative di- 
electric constant £i- are position-dependent material pa- 
rameters in a general three-dimensional nanostructure 
that may be composed of several different materials. The 
total charge density is composed of the electron density 
n (x) and the ionized donor density (x) , 



p(x)--n(x) + Ar+(x) 



(3) 



The latter results from the donor density A'd (x), the 
degeneracy ^d, and the donor energy iJo according to 
the Thomas-Fermi approximation. 



AfD(x) 



1 + qj3e(^F--^^D)/'^B 



(4) 
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The electron density is calculated quantum mechanically 
by occupying the electronic states according to the Fermi- 
Dirac statistics 

n(x)=^|(x|a,)|V(^^^^^). (5) 

Here, T denotes the temperature, fee is Boltzmann's con- 
stant, E-p denotes the Fermi level, and / is the Fermi 
distribution function. 



tion matrix elements 
{a,Y\Vi2\ajZ) 

(11) 

where IctiK) is the product state of particle 1 state \ai), 
and particle 2 state \Y). 



B. Electronic structure of particle 2 

For the second particle, we will focus on an electron 
in a closed double quantum dot that can tunnel between 
these dots. The two lowest lying quantum states can be 
described by an effective Hamiltonian, that takes into ac- 
count the energy splitting A between the ground states 
of the isolated quantum dots and the tunneling coupling t 
between the two quantum dots. This two- by- two Hamil- 
tonian is given by 



H2\Y) = 



A 



^2 \Y)=Ey\Y). 



(6) 



Typical values for A and t for lateral semiconductor 
quantum dots are of the order of 10 fieVM We use this 
model Hamiltonian for the subsystem of the second par- 
ticle and represent the charge distribution of the elec- 
tron in the ground state of each (isolated) quantum dot 
by a point charge centered at Xq and Xi, respectively. 
The use of a model Hamiltonian rather than a realistic 
three-dimensional device Hamiltonian is not a principle 
limitation of the present method but is adequate for the 
concrete device geometry that we will study in detail in 
this paper. The eigenstates of H2 are linear combina- 
tions of the ground states |0) and |1) of the two isolated 
quantum dots, 



\B)=hBo\0) + hBi\l), 
\A) =hAo\0) + hAi\l), 



(7) 

(8) 



with real valued coefhcients hyj (Y — A,B, J ~ 0,1). 
Here, \B) is the bonding eigenstate and \A) the anti- 
bonding eigenstate. The corresponding eigenenergies are 



EA^^Vt^ + A^- 



(9) 
(10) 



For vanishing t, the eigenstates \B) and \A) reduce to |0) 
and 11). 



C. Interaction matrix elements 

The Coulomb interaction V12 between particle 1 and 
particle 2 yields the following expression for the interac- 



D. Ballistic transport through a system of two 
interacting particles 

In this second step, we develop a Green's function 
method by extending the CBR method^i^i to the case of 
an open device that describes a system of two interacting 
particles. We stress that we will consider the two parti- 
cles as distinguishable. This is not an approximation for 
the proposed two-qubit devices but comes directly from 
the system geometry. Indeed, we choose the QWRs and 
the double dot to be well separated from each other so 
that no significant tunneling between the two structures 
can occur. In other words, the two wave functions - 
namely the bound state in the double dot and the scat- 
tering state in the QWRs - never overlap. This leads to 
the distinguishability of the two particles based on their 
spatial localization. In addition, we assume that the in- 
teraction is negligible outside the device, due to screening 
in the contacts and/or due to barriers. In the following, 
the term device denotes a finite three-dimensional region 
that is coupled to reservoirs by an arbitrary number of 
leads. The device may be under applied bias and contain 
some spatially varying potential. The total (two-particle) 
Hamiltonian of the system, including the device and the 
leads, can be written in symbolic matrix form 



Htot = 



■•• 



\wl ■■■ 







W-, 



Nl 



Wn, 



(12) 



where represent the Hamiltonian of lead A, the 
Hamiltonian H'^ corresponds to the device region, and 
W\ is the coupling between the device and this lead 
(A — 1, . . . , Nl). H° is composed of the single-particle 
Hamiltonians Hi and H2, corresponding to particle 1 and 
2, respectively, and the interaction term V12, 



= Hi+H2 + Vi 



12- 



(13) 



The leads (acting as reservoirs) are semi-infinite and 
therefore, the total Hamiltonian i?tot is infinite- 
dimensional. This infinite-dimensional Hamiltonian can 
be reduced to a non-Hermitian finite-dimensional Hamil- 
tonian H = -t- E that describes the open device 
exactly^. In this formulation, the influence of the 
leads is included through a finite-dimensional operator 
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E = Ei + -- - + Sjv£,- This is the sum of the complex con- 
tact self-energies E^ that are nonzero only in the contact 
regions where the lead A adjoins to the device. The self- 
energies couple only particle 1 to the leads since particle 2 
is assumed to be bound. The Hcrmitian Hamiltonan iJ" 
represents the decoupled device, i.e. the device with no 
coupling to the leads. Hi as well as H2 are Hamiltonians 
of the closed device. In the ballistic case, all observables 
of interest such as transmission functions and the current 
can be calculated from the retarded Green's function 
of the open device. This is defined by 



{E - H)-^ = {E-Hi-H2- V12 - E)-i . (14) 



The retarded Green's function can be obtained from the 
Green's function of the decoupled device by the Dyson 
equation 



(15) 



which in turn can be evaluated from its spectral repre- 
sentation 



G" 



\n) {n\ 



(16) 



(17) 



The direct evaluation of G^ according to Eq. (jlSp re- 
quires the inversion of a huge matrix that is proportional 
to the number of grid points Nu of the device. By con- 
trast, the GBR method allows one to drastically reduce 
this effort by utilizing the following exact properties of 
G^ that remain valid in the present many-particle case. 
We decompose G^ into four blocks, 



G'' = 



^CD 



^DC 



(18) 



where G§ is a matrix proportional to the number Nc of 
contact grid points and is called the contact block. Note 
that Nc ^ Njj. For the contact block G^, the following 
equation holds: 



G« = 



If 



1 r'O 



(19) 



where l^, Gq, and E^ are the corresponding contact 
blocks of the unity matrix, G*^, and E, respectively. The 
key point is that this is a linear matrix equation of the or- 
der of Nc, since E is nonzero only in the contact region. 
As a consequence, this contact block is also sufficient 
to determine the transmission functions through the de- 
vice. The ballistic current can then be calculated from 
the transmission functions according to the Landauer- 
Biittiker formalism. In the following, we represent G^ 
in the mixed basis of position eigenstates of particle 1 
and energy eigenstates of particle 2, 

G«(x„x„r,Z) = Yl {^^\a^){a^Y\G^\a,Z){a,\^,). 

(20) 



In this form, G^ represents the probability amplitude 
for the propagation of particle 1 from position x^- to po- 
sition Xi, accompanied by a transition of particle 2 from 
eigenstate \Z) to For A 7^ A', the expression 



T^;f{E) = Trir^(Gg)r^' G 



. (e^ ~ E^t) . 



(21) 
(22) 



therefore yields the probability for the transmission of 
particle 1 from lead A' to lead A under the condition, 
that initially particle 2 is in state \Z) and ends up in state 
\Y). Note that the trace in Eq. (|2ip is taken only with 
respect to the position variables of particle 1 . If the final 
state of particle 2 is not measured, the corresponding 
transmission probability is obtained by summing over all 
final states of particle 2 



AA' 



Y,Z 
AA' ■ 



(23) 



In our concrete calculations (see Sec. |V]), we found it 
most efficient to first determine the eigenstates of the 
single-particle Hamiltonians Hi and H2 and then diag- 
onalize the interaction V12 in the product basis of these 
single-particle eigenstates. This procedure allows us to 
take advantage of another built-in efficiency of the GBR 
method, namely the fact that only energetically low lying 
eigenstates of the closed device couple to the incoming 
and outgoing lead statesi^i^ Gonsequently, it suffices to 
take into account only a reduced set of eigenstates of the 
single-particle Hamiltonians Hi and H2 without affecting 
the transmission results noticeably. 

Another quantity that can be readily obtained from 
the retarded Green's function G^ is the charge density 
of the current carrying states of particle 1. This requires 
one to know not only the contact block G^ but also the 
submatrix G^c- However, the calculation of the subma- 
trix G^c only requires one to evaluate the inverse of the 
small matrix Ic — EpG^, 



^R. _ r<0 
'-'DC — '-'DC 



1(7 — Gr<I^ 



c^cj 



(24) 



where the matrix G^^ corresponds to G|^p in the open 
device.-2i According to Eq. (j20p . the latter Green's func- 
tion reads in the mixed basis 

Ggc(x.,x,,r,Z) = {^,\a.^{a.,Y\GicW,Z){a,\^.j). 

(25) 

The total charge density of the current carrying states of 
particle 1 is consequently given by 



(26) 
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where the contribution 

(27) 

X r^(x,,x,) (Ggc(x.,x„r,Z))^ (28) 

results from the projection of the total charge density 
onto the eigenstate \Z) of the second particle. Here, 
denotes the Fermi level in lead A. 



III. NUMERICAL DETAILS 
A. Electronic structure of particle 1 

Equations ([!]) and ([2]) are coupled partial differential 
equations in position space that we discretize and map 
onto a nonuniform tensor product grid. For the dis- 
cretization of the all-electric Mach-Zehnder interferom- 
eter, a total of 6 X 10^ grid points have been used. We 
employ box integration finite differences in order to en- 
sure flux conservation across boundaries with different 
material parameters. For the solution of the nonlinear 
system that results from the discretization of the nonlin- 
ear Poisson equation Newton's method with inexact 
line search is invoked. This remaps the problem into a 
sequence of linear solutions. We use the Dupont-Kendall- 
Rachford preconditioned conjugate gradient method for 
the solution of the resulting linear systems. -■^'^ The dis- 
cretization of the Schrodinger equation (jT]) results in a 
large matrix eigenvalue system. However, since the oc- 
cupation of electron states falls off exponentially with in- 
creasing energy distance from the Fermi level, only about 
250 quantum states suffice for a converged calculation of 
the electron density. This allows for the use of the itera- 
tive Arnoldi method that is implemented in the published 
ARPACK librariesJS We found that the calculation of 
eigenvalues and eigenstates can be accelerated by an or- 
der of magnitude by invoking a Chebyshev polynomial 
based spectral transformation that provides an efficient 
preconditioning of the linear system4i Finally, in order 
to solve the coupled Poisson-Schrodinger system, we em- 
ploy an approximate quantum charge density inside of 
Poisson's equation in order to estimate the dependence 
of the density on the potential through Schrodinger's 
equation. This corrector-predictor scheme reduces the 
number of required diagonalizations and accelerates the 
convergence of the coupled system significantly4^ 

Equation Q requires the specification of boundary 
conditions in order to obtain unique solutions. For the 
potential in the Poisson equation, we generally use von 
Neumann boundary conditions. Dirichlet boundary con- 
ditions are employed at Schottky contacts and at the 
GaAs/air interfaces. In the former case, the electrostatic 



potential at the boundary is set equal to the Schottky 
barrier height, i.e. the difference between the conduction 
band edge and the quasi- Fermi level. At the GaAs/air 
interface, on the other hand, the electrostatic potential 
is determined by the experimentally known Fermi level 
pinning. The Schrodinger equation for the open quan- 
tum system is solved in terms of the CBR method that 
has been described in detail elsewherej2i For the calcu- 
lation of the transmission probabilities, we have taken 
into account the lowest 250 eigenstates for the determina- 
tion of the retarded Green's function of the closed device 
and the lowest 5 propagating modes in each of the four 
leads. The self-consistent solution of Eqs. ([T]) and ([2]), for 
given split gates voltages, takes approximately 2 days on 
state-of-the-art PC hardware. Converged self-consistent 
results require typically 12 iterations in the two sets of 
equations. 



B. Entangled Mach-Zehnder double quantum dot 
device 

Our novel method for the calculation of the ballistic 
transport properties through a system of two interacting 
particles requires the diagonalization of the two-particle 
Hamiltonian including the interaction. This has to be 
done once for every bias point in the 1-V characteristics. 
The diagonalization is performed in the set of product 
states formed by subsets of the single-particle eigenstates 
of the two non-interacting particles. We found that a 
subset of 100 lowest eigenstates of particle 1 is sufficient 
to obtain convergent transmission functions in the meV 
energy range around the Fermi level that is relevant for 
this device. 



IV. ANALYTICAL MODEL OF THE 
TWO-QUBIT DEVICE 

Before presenting and discussing the numerical results, 
we develop a simple analytical model of our two-qubit 
Mach-Zehnder double-dot system in the present section. 
The device is shown schematically in Fig.[T] The purpose 
of this procedure is to be able to better grasp the physics 
behind the generation of entanglement before we present 
the numerical results of the fully three-dimensional self- 
consistent calculations in the next section. 

The good agreement between the two approaches will 
be a further confirmation of the correct functioning of 
the device: The results that we obtain for the entangled 
Mach-Zehnder double quantum dot device are indeed the 
result of a two-qubit quantum operation. Obviously, the 
results of the two methods will only be in qualitative 
agreement with one another, since the analytical model 
developed in the this section cannot account for the full 
complexity of the realistic device. 

We describe the device consisting of a Mach-Zehnder 
interferometer coupled to a double quantum dot, by two 
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FIG. 1: (Color online) Schematic view of the proposed quan- 
tum transport device. The device is realized by two stacked 
GaAs/AlGaAs 2DEGs. The top 2DEG (green) is depleted by 
external gates to form a double-QWR based Mach-Zehnder 
interferometer. The bottom 2DEG (blue) contains two elec- 
trostatically defined, coupled quantum dots. The larger one is 
centered beneath and in between the QWRs. There is a bias 
voltage VsD applied between the source (5*) and the three 
drain contacts (Di — D3) that leads to currents (Ji, J2) that 
are labeled accordingly. There are additional gates acting on 
the device that are shown in the following Fig. [2] For sake of 
clarity, the figure is not drawn to scale. 



coupled qubits, W and D, respectively. The basis states 
|0)^ and |1)^ of qubit W are defined by an electron 
that propagates through either of the two QWRs. The 
basis states |0)^ and |1)^ of qubit D, on the other hand, 
are defined by the ground states of the isolated quan- 
tum dots. Quantum gates can be described by unitary 
operators that map an initial state |$in) onto a corre- 
sponding final state |$fi). The unitary operator T that 
represents our two-qubit quantum transport device can 
be decomposed into a sequence of unitary operators of 
less complex single- and two-qubit gates as follows^: 

f = RpV-R. (29) 

The operators R and P represent single-qubit gates that 
act solely on qubit W. In particular, R represents a 7r/2 
rotation gate and describes a coupling window, whereas 
P represents a phase gate that adds a phase difference 
to the two components of the qubit. The latter can be 
realized through an additional split gate located on one of 
the QWRs, as explained in Sec.|Vl The sequence R P R 
therefore corresponds to a Mach-Zchnder interferometer. 
The operator V represents a two-qubit gate and describes 
the coupling of the two qubits D and W due to their 
interaction. This decomposition allows us to specify the 
action of the single-qubit gates R and P directly in terms 
of unitary operators. 

We proceed by describing the entangled Mach-Zchnder 
double dot device as a stationary quantum-mechanical 
scattering problem. To this end, we represent the two 
electrons in the pair of QWRs and the double quantum 
dot by the following 8 basis states for given total energy 
E and a given position a;, 

'fp±k,,aj (x) = exp {±ikjx) |ct)^ | J)° . {cr,J ^ 0, 1) 

(30) 

These states describe the two electrons as follows: The 
first electron propagates along the wire axis x in either 



of the two states |0) and |1) , corresponding to the 
two QWRs. The electron momentum at position x is de- 
noted by kj. The second electron occupies one of the two 
quantum dots, represented by states jO)'^ and In 
this representation, the operators in Eq. (|29p correspond 
to 8 X 8 transfer matrices, 

T = Tf^-Tp-Ty-T^. (31) 

This will allow us to finally determine the transmission 
function through the device. As a first step^ we derive the 
transfer matrix Ty of the two-qubit gate V by solving a 
two-particle Schrodinger equation. The non-interacting 
two-particle Hamiltonian that describes the dynamics of 
the system is given by 

Ho = ^+H2. (32) 

Here, H2 is the model Hamiltonian that describes a pair 
of tunneling coupled quantum dots in terms of the bare 
splitting A and the tunneling coupling t in accordance 
with Eq. For the Coulomb interaction between the 
two electrons, we make the following ansatz 

V (x) = [/ |1)^ \lf (1|^ (1|° (x) 9{L-x). (33) 

This describes a localized interaction with interaction 
strength U that is nonzero only for the two-particle state 
|1)^ |1)^ but zero for all other states. In addition, the 
interaction vanishes if the position of the QWR electron 
lies outside the spatial interval from x = to x = L. The 
resulting quantum mechanical scattering problem 

H = Ho + V{x), (34) 

can be solved analytically for the 8x8 transfer matrix 
Ty{E). Since V (x) only acts on the state |1)^ |1)°, Ty 
has the following form, 

Ty = C'q^' , (35) 

where r is a 4 x 4 submatrix that connects the states 
a — 1, i.e. states that contain |1)^. In order to deter- 
mine this submatrix, we split the a;-axis into 3 regions 
x < 0, < X < L, and x > L. Within each region, the 
Hamiltonian H can be diagonalized. Requiring continu- 
ity and differentiability, we obtain the following explicit 
expression for the submatrix r. In this expression, we 
have transformed the basis by using the bonding and an- 
tibonding dot states \A)^ and \B)^ from Eqs. O, © 
instead of the states |0)° and |1)^. In this basis, we ob- 
tain r = nj^^ iij^Kq^ iiQ where the submatrices fi, k read 
for xi = Q, X2 = L 
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where i = 1,2. The 2x2 matrix 1/^^ (k) is defined as 
foUows: 



_ f exp (ikxi) exp {—ikxi) 
ikexp{ikxi) —ikexp{—ikxi 



(38) 



The coefficients ay, and the wave vectors fcy, qy 
{Y = A, B) are given by [t, U > {)) 



Uy 



tu 



{AEy {Ey +e) + AUf + t^U^ 



AEy {Ey + £) + MJ 
{AEy {Ey +e) + AUf + t^U^ 



(39) 



kY = -\/2m {E~Ey), 



.Y^'hrr^fE-luTs 



(41) 



(42) 



where e 



iyt2 + ([/ + A) . Note that the electron 
wave vectors fcy arise from the region a; < and from 
X > L, whereas qy arises from the central region. 

The next step is to determine the transfer matrices 
and Tp from the following representations of ideal 
rotation and phase gates, 



R = 



cos{e/2) 



{em\ 



sin (0/2) cos(6'/2) / ' 



P = 







-10/2 



= vr/2, (43) 



(44) 



where </> is the phase gate angle. We assume that these 
gates do not introduce back-scattering. In this case, we 
can apply R and P separately for both directions of prop- 
agation. This leads to unitary matrices Tp,Tp with the 
following nonvanishing elements 

(T^),, , = cos {9/2) , (T^)m+4 = (-1)' sin(0/2) , 
(r^),+4,. = (-l)*+^sin(0/2), (j = l,4;j = l,8) (45) 



{Tp)^^^ = exp 



^ ' 2 



{j = l/' 



(46) 



This completes the calculation of the total transfer ma- 
trix T. The final step is to determine the transmission 
function which requires us to set an initial condition for 
the scattering problem. We consider the situation where 
the propagating electron enters the pair of QWRs from 
the left in state |0)^. We further assume that the elec- 
tron in the double quantum dot lies initially in the bind- 
ing state \B)^ . The eight final state coefficients rep- 
resent the transmission amplitudes into the four states 



|cr)^|y)° (cr = 0,1, F = A,B) for each of the two 
asymptotic propagation directions (to the right or to the 
left). We denote these final state coefficients by E^y 
and i?^ y , respectively. By construction, the two rota- 
tion gates and the phase gate do not introduce any back- 
scattering. The amount of back-scattering due to the 
interaction V {x) depends on the ratio of the interaction 
strength U and the length L of the interaction region 
and can be made arbitrarily small. Thus, the transmis- 
sion amplitudes E"^ y are negligible. In addition, we only 
need the transmission probabilities Ti and T2 that de- 
scribe the total probability of an electron to arrive at 
|0)^ or |1)^, irrespective of the state of the double dot 



(40) electron, 



Ti — \Eo^a\ + 1^0, b\ 
T2 = I^^I.aT + I^IbI 



(47) 
(48) 



The resulting transmission probabilities Ti and T2 will 
be presented in the following section for different values 
of the tunneling coupling. They will be compared with 
the numerical results obtained with the realistic model 
described in Sec. HIl Indeed, the above comparison will 
show that the two models are able to reproduce the ba- 
sic physics of the device and that the device is able to 
produce an entangled state of the two qubits. 



V. RESULTS AND DISCUSSION 

The method presented in the Sec. HIl allows us to calcu- 
late the ballistic current through entangled two-particle 
systems. Concretely, we first present results for a double- 
QWR based Mach-Zehnder interferometer and show that 
this interferometer can be employed as a fully control- 
lable single-qubit gate. Based on this system, we then 
propose a novel two-qubit quantum transport device that 
allows the controlled generation and measurement of 
entanglement between the QWR qubit and the double 
quantum dot qubit. The measurement involves only the 
DC I-V characteristics, no higher order current correla- 
tions or magnetic fields are required to detect entangle- 
ment. We would like to emphasize that the subject of 
this paper is to show that these devices are well suited to 
prepare arbitrary single-qubit and entangled two-qubit 
states. Actual quantum computing applications, how- 
ever, require additionally the performance of these oper- 
ations in a time-resolved fashion with single electrons. 

A schematic view of the proposed two-qubit device 
is shown in Fig. [TJ It is based on a GaAs/AlGaAs 
heterostructure that consists of two vertically stacked 
2DEGs. The 2DEGs are locally depleted by the appli- 
cation of negative voltages to appropriate metal top and 
bottom gates, that are not shown in this figure. For suit- 
able gate voltages, the top gates create a pair of parallel 
QWRs in the upper 2DEG that are connected to each 
other by two coupling windows. For specific energies, 
these coupling windows act as beam-splitters and the 
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QWR network behaves as an electronic Mach-Zehnder 
interferometer. In contrast to electronic Mach-Zehnder 
interferometers based on quantum Hall edge channels,— 
no magnetic fields are employed. We therefore use the 
term all- electric Mach-Zehnder interferometer for this 
kind of device. We note that the relative phase shift 
between the electronic wave function in the two QWRs 
is controlled electrostatically by phase gates. These are 
additional top gates, but may be combined with the gates 
that are used to define the QWRs in the first place. The 
lower 2DEG contains the two electrostatically defined 
tunneling-coupled single-electron quantum dots. An elec- 
tron in the upper 2DEG couples to the electron in the 
double quantum dot in the lower 2DEG by means of the 
Coulomb interaction. We have designed the entire de- 
vice so that the consequences of entanglement become 
markedly evident in the observables, and we refer to it as 
Entangled Mach-Zehnder Double Quantum Dot Device. 

The Mach-Zehnder interferometer is operated under 
ballistic transport conditions. By applying a small DC 
bias voltage Vsd between the upper left source S and the 
remaining three drain contacts Di, D2, -D3, the device can 
be used in such a way that current flows predominantly 
from S to Di (denoted by Ji) or to D2 (denoted by J2), 
as indicated in Fig. [TJ 



A. Results: All-electric Mach-Zehnder 
interferometer (single qubit device) 

First, we turn to the I-V characteristics of the single- 
electron device, i.e. the all-electric Mach-Zehnder inter- 
ferometer alone. The calculations have been based on the 
GaAs/AlGaAs heterostructure that is depicted schemat- 
ically in Fig. [21 It consists of a 5 nm thick GaAs cap 
layer [see vertical cross section in Fig. [2] (a)], followed by 
a 45 nm barrier of Alo.37Gao.63As, and a 10 nm GaAs 
2DEG layer. Beneath, there lies a 2 /im Alo.37Gao.63As 
substrate. A silicon S doping layer with a concentration 
of 2.5 X 10^^ cm~^ is located 25 nm below the surface. 
The electron sheet density in the 2DEG layer has been 
calculated to be 2.4 x 10"'^^ cm~^ for the ungated sample 
at a temperature of 4 K. 

Figure [2] (b) shows a top view of the structure includ- 
ing the gates. The source and drain contacts (cf. Fig. [T]) 
are not shown in this figure. The device is 800 nm wide 
and 1345 nm long. Here we distinguish three types of 
gates: side gates (black), mid gates (gray), and phase 
gates (hatched). Each of the two side gates are, respec- 
tively, 280 nm wide, 1345 nm long, and are biased at 
VsG (S) -0.245 V with respect to the source. The 
three mid gates are 40 nm wide and 200 nm, 800 nm, 
and 200 nm long, respectively. We apply a gate voltage 
of V"sG (M) = -0.660 V to the mid gates. They define two 
QWRs (white) with a nominal width of 100 nm as well 
as two coupling windows, each with a nominal length of 
72.5 nm in the direction parallel to the wires. All of these 
structural parameters have been chosen to guarantee op- 




FIG. 2: (a) Vertical cross section of the GaAs/AlGaAs het- 
erostructure with a single 10 nm wide GaAs quantum well 
that is used for the Mach-Zehnder interferometer. GaAs is 
shown in black, Alo.37Gao.63 As is shown in gray. The GaAs 
cap layer is 5 nm thick. A silicon 5 doping layer with a concen- 
tration of 2.5 X lO'^^ cm~^ is located 25 nm below the surface 
(dashed line), (b) Top view of the Mach-Zehnder interferom- 
eter that depicts the side gates (black) with voltage Vsg(S) 
and the mid gates (gray) with voltage Vsg(M) that define the 
two QWRs and the coupling windows between them. On top 
of the wires, there are two additional phase gates (hatched) 
at a bias of ±Vsg(P)- 



timal device operation, while allowing for its fabrication 
with current technologies. 

The subband spacing of the two lowest subbands in 
each of these QWRs amounts to 3.1 meV. It is impor- 
tant for the interferometer that the Fermi wave length is 
close to the length of the coupling windows. Indeed, our 
calculations yield a Fermi wave length of Ap = 77 nm 
in the lowest subband. In order to control the relative 
phase of the electron wave function in the two QWRs, a 
small gate voltage -I-Vsg (P) and — Vsg (P) is applied to 
the phase gates, respectively. The two phase gates are 
centered on top of the two QWRs and are each 720 nm 
long and 95 nm wide. A separating layer of resist could 
be used to reliably insulate the phase gates from the other 
gates. Figure [3] shows the equilibrium charge density (up- 
per panel) and the corresponding potential (lower panel) 
in the 2DEG layer of the Mach-Zehnder interferometer. 
The black and the white framed rectangles indicate the 
position of the metal top gates. The upper panel shows 
the 2DEG to be fully depleted underneath the gates. In 
addition, it is evident from the figure that the QWRs and 
the coupling windows are formed indeed. The wires are 
strictly one-dimensional in the sense that only the low- 
est subband contributes to the density. The lower panel 
shows the potential barrier that separates the two QWRs. 
The energy scale is chosen so that the Fermi level lies at 
zero meV. 

Figure[4](a) shows the currents Ji (solid), J2 (dashed), 
as well as the back-scattered current J3 (dotted) that 
flow from the source S to Di, D2, and D3, respectively. 
The currents have been calculated by assuming a DC 
bias voltage of 50 and a temperature of 30 mK. The 
figure shows the dependence of the currents on the volt- 
ages -|-VsG (P) and —Vsg (P) applied to the two phase 
gates, respectively. We first note that the back-current 
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FIG. 3: (Color online) Top views of the equilibrium electron 
density (upper panel) and electrostatic potential (lower panel) 
within the upper 2DEG. The black and the white frames in- 
dicate the position of the metal side and mid gates. 



J3 is smaller than 0.5 nA for the entire relevant range 
of gate voltages and therefore negligible. In contrast to 
J3, the currents Ji and J2 strongly depend on the gate 
voltage and oscillate between and 4 nA for more than 
1.5 oscillation periods. The currents Ji and J2 are seen 
to be phase-shifted relative to one another by 180 de- 
grees which confirms that the device actually behaves as 
a Mach-Zehnder interferometer. Note that this Mach- 
Zehnder interference pattern becomes damped out once 
the absolute gate voltage \Vsg (P) I exceeds 15 mV be- 
cause this increases the amount of back-scattering. Back- 
scattering can be suppressed, however, by increasing the 
length of the device, because the phase shift depends 
linearly on the length of the phase gate, whereas back- 
scattering is independent of the length but increases with 
the magnitude of |Vsg (P) |. 

For comparison, Fig. |4] (b) depicts the calculated cur- 
rents Ji , J2 , and J3 for the same device but for a nominal 
length of the coupling windows that has been reduced by 
5 nm. The difference in results shows the sensitivity of 
the interference pattern to small changes in geometry, 
and implies very high demands on fabrication precision. 

Figure [5] shows the charge densities of the station- 
ary current carrying states of the Mach-Zehnder inter- 
ferometer for selected phase gate voltages. The three 
Figs. [5] (a), (b), and (c) correspond to gate voltages 
of Vsg(P) = mV, —3.6 mV, and —7.5 mV, respec- 
tively. The three graphs show that the quantum me- 
chanical charge densities change predominantly near the 
drain contacts on the right hand side of the device. These 
changes reflect the redistribution of the total current be- 
tween Ji and J2. In Fig. [5] (a), the charge density is zero 
near the upper right terminal but large near the lower 
right terminal. This agrees with Fig. [4] (a) that shows 
that the current Ji is minimal ( J2 maximal) for the cor- 
responding gate voltage Vsg(P) = mV. In Fig. [5] (b). 
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FIG. 4: (a) Stationary currents Ji (solid), J2 (dashed), and 
J3 (dotted) through the Mach-Zehnder interferometer as a 
function of the phase gate voltage Vsg (P) for a nominal length 
of the coupling windows of 72.5 nm. The source-drain voltage 
VsD is set to 50 /iV. (b) Same as (a), but for a coupling 
window length of 67.5 nm. 



on the other hand, the charge density is seen to be al- 
most equal near both right drain contacts, again in ac- 
cord with Fig.H (a) for Vsg(P) = -3.6 mV. Finally, the 
charge density is large near the upper right terminal but 
zero at the lower right terminal in Fig. [5] (c), again in 
correspondence with Fig. 2] (a). 



All of these effects are a consequence of the self- 
interference of the electron wave function. Depending 
on the relative phase between the two QWRs, the in- 
terference at the second beam splitter leads to partial 
extinction either near the upper or the lower right termi- 
nal. 
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FIG. 5: (Color online) Charge densities associated with the 
current carrying scattering states for a phase gate gate voltage 
VsciP) of (a) meV, (b) -3.6 meV, (c) -7.5 meV. In all cases, 
the source-drain voltage Vsd is set to 50 fiV. 
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FIG. 6; Schematic representation of the Mach-Zehnder inter- 
ferometer. We have indicated the three quantum-operation 
regions, consisting of the two coupling windows R and one 
phase gate P, and four other segments where the qubit state 
is well defined (see Table |I|. 



B. Discussion: Switching characteristics of the 
all-electric Mach-Zehnder interferometer 

The DC transfer characteristics shown in Fig. |4] (a) ex- 
hibits multiple pronounced maxima and minima that can 
be attributed to rotations of a qubit on the Bloch sphere. 
Before entering into the discussion of the results, let us 
give a rigorous definition of the Mach-Zehnder interfer- 
ometer qubit state. To this aim, we depict the device 
schematically in Fig. [5] where we show the regions of the 
two coupling windows and the phase gate that perform 
the R and P transformations of Eqs. (|43| and (|44p. re- 
spectively. We define the basis states |0) and |1) of the 
QWR qubit by an electron scattering state localized in 
the upper and lower QWR, respectively. This definition 
holds only in the regions where the two channels are well 
separated, as in the segments I, II, III, and IV of Fig. [S] 
By contrast, the qubit state is not well defined within 
the window regions R, where the barrier between the 
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TABLE I; Quantum state of the Mach-Zehnder qubit in the 
four device segments indicated in Fig. [6] (columns), corre- 
sponding to the three situations shown in Fig. [5] (rows), re- 
spectively. 
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FIG. 7: Contour plot of the transmission function Ti between 
source S and drain Di as a function of kinetic energy in meV 
and as a function of the coupling window length L in nm. The 
total width of the structure equals W — 220 nm, whereas 
the width of the inner coupling windows amounts to W' = 
100 nm. The darker the color, the smaller the value of Ti. 
The dark lines therefore indicate minima in Ti. The arrows 
point to crossings of resonances where Ti vanishes. 



QWRs is very small. As an example, the approximate 
qubit states for the four regions indicated in Fig. [6] are 
included in Table HI and specified for each of the three 
cases depicted in Fig. O Indeed, in the calculation of 
Fig. [5l the split gate voltages have been chosen in order 
to obtain a qubit rotation angle 9 = n/2 and phase angles 
(j) = 0,7r/2,7r in (a), (b), (c), respectively [cf. Eqs. ([^5)1 
and ((44| ]. 

Since it is unlikely that one can fabricate the device 
discussed in Fig. [2] precisely, it is imperative to under- 
stand the dependence of the Mach-Zehnder interferences 
on the various device geometry parameters such as width 
and length of the QWRs and the coupling windows not 
only quantitatively but also qualitatively. A physically 
transparent picture can be obtained by studying a model 
H shaped structure that consists of a single coupling win- 
dow where the depletion potentials are hard-wall poten- 
tials as depicted in the inset of Fig. \7\ 

The coupling between the two QWRs causes the trans- 
mission to become resonant for specific energies that de- 
pend on the length L and width W of the coupling win- 
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FIG. 8: Transmission functions Ti (solid), T2 (dashed), andTs 
(dotted), as a function of the kinetic energy of the electron in 
meV. The length L of the coupling window equal has been set 
to 82 nm. The other symbols have been defined in Figure [T] 
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FIG. 9: Currents Ji and J2 in nA as a function of the phase 
gate voltage Vsg(P) in mV for two different quantum dot 
tunneling couplings t in /leV. The inset is a simplified version 
of Figure [T] 



dow. These energies are approximately equal to the en- 
ergies of the bound states in a rectangular box with the 
dimensions of the coupling window, i.e. 

Ti 711 

E{n,m)oc— + —. (49) 

Beam splitting, i.e. a crossover in the transmission from 
channel S-Di to S-D2, requires that the correspond- 
ing transmission functions obey the relation Ti{Eq) w 
T2{Eq) « 0.5 for some energy Eq. In turn, this condition 
requires that there are 2 sets of values {n, m) that corre- 
spond to the same energy i.e. E(ni, mi) — E{n2, m2). 
Thus, the crucial geometric requirement for quantum me- 
chanical switching is to achieve an appropriate ratio of 
length L and width W that obeys this condition. In or- 
der to obtain a pronounced resonance, the values of n, m 
should be of the order of unity additionally, since the 
resonances lie densely in energy for large integers. Fig- 
ure [7] shows a contour graph of the transmission proba- 
bility Ti (from source S to drain Di) as a function of the 
coupling window length L and the kinetic energy of the 
electron. The arrows in this figure point to suitable val- 
ues of the length L where the energies of two eigenstates 
cross. To further illustrate this crossover in the transmis- 
sions. Fig. [5] shows a cross-section through Fig. [7] for the 
length L = 82 nm that is marked by the leftmost white 
arrow. This crossing corresponds to the switching from 
channel S-Di to S-D2 and occurs near the degeneracy of 
the energy states £'(1,3) and E{2, 1). 

C. Results: Entangled Mach-Zehnder double 
quantum dot device (2-qubit device) 

In the following, we present the results for the entan- 
gled Mach-Zehnder double quantum dot device that has 
been schematically depicted in Fig. [T] For these calcula- 
tions, we have used a somewhat different set of geometry 



parameters for the Mach-Zchnder interferometer in or- 
der to optimize the entanglement. Since the calculations 
of the entangled system are very demanding, we have 
defined the Mach-Zehnder interferometer by a hard-wall 
potential within a 10 nm thick slab of GaAs rather than 
performing a charge self-consistent calculation including 
all metal gates. We have taken the two QWRs to be 
55 nm wide, 1000 nm long, and the lateral distance be- 
tween them has been set to 20 nm. The coupling windows 
now have a length of 85 nm. The resulting energy spacing 
of the two lowest subbands amounts to 5.3 meV which 
is nearly twice as large as in the Mach-Zehnder device of 
Sec.|V]A. The Fermi level Ep has been set to 1.6 meV in 
the lowest subband. This Fermi level causes the two cou- 
pling windows to act as almost perfect beam-splitters, 
with channel transmissions Ti{Ep) « T2{Ep) « 0.5. 
For the calculation of the interaction between the Mach- 
Zehnder interferometer and the double quantum dot, we 
assume a vertical distance of 80 nm between them. The 
larger of the two quantum dots with a lower ground state 
energy is located exactly underneath the center of the 
Mach-Zehnder interferometer whereas the smaller one 
lies 60 nm away in the direction perpendicular to the 
QWRs. The double quantum dot system is modeled by 
the Hamiltonian given in Eq. ([6]). The positions of the 
Mach-Zehnder interferometer and the double quantum 
dot relative to each other have a large influence on the 
operation characteristics of the device. The present con- 
flguration exploits the idea that the interaction between 
the electron in the QWRs and the electron in the larger 
quantum dot is identical for both wires. The electron in 
the smaller quantum dot, on the other hand, lies closer 
to the QWR between gates D2 , D3 and interacts mostly 
with the electron in that wire. Consequently, a phase dif- 
ference in the electron wave function for the two QWRs 
is exclusively induced by the electron in the smaller dot. 
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FIG. 10: Visibility as defined in the text, corresponding to the 
interference pattern of Fig. |9l as a function of the tunneling 
coupling in /^eV. 



At first, we study the transfer characteristics of the 
device for a situation where the tunneling between the 
quantum dots is inhibited {t = 0). The ground states 
of the larger and the smaller quantum dot are assumed 
to differ by A = 10 /zeV. The resulting current in the 
QWRs is shown in Fig. [5] by the dashed curves. We ob- 
tain the typical interference pattern that we have also 
found in Sec.|V]A for the all-electric Mach-Zehnder in- 
terferometer. The asymmetry with respect to the sign 
of the gate voltage is caused by the asymmetric position 
of the phase gate along only one of the two QWRs in 
combination with the additional repulsive potential due 
to the electron in the double quantum dot. 

We now allow the electron in the double quantum dot 
to tunnel between the two dots {t = 10 iieV). The cur- 
rents Ji and J2 that correspond to this situation are 
shown in Fig. [5] by the solid curves. We still obtain the 
Mach-Zehnder interference pattern, but the visibility is 
markedly reduced. Let us define the visibility by 

(J2- Ji)/(J2 + Ji)|^.3^(P)=o' (50) 

i.e. for zero phase gate voltage Vsg(P)- Figure [TO] shows 
the visibility as a function of the tunneling coupling. The 
visibility is almost 1 for vanishing tunneling coupling and 
decreases monotonously with increasing tunneling. The 
visibility is therefore uniquely related to the magnitude 
of the tunneling coupling. The above effect can be ex- 
plained as follows. If the tunneling probability is high, 
the double-dot electron can easily change its state. In 
this regime, it is also easy for the double-dot to measure 
the path of the interferometer particle, thus causing de- 
coherence and destroying the interference. For low tun- 
neling probability, on the other hand, the electron that 
sits initially in the large dot is insensitive to the Coulomb 
interaction with the interferometer. When the tunneling 
is completely suppressed, the state of the dot is fixed 
and the interferometer electron is unable to change the 
dot state. While the electrons are still interacting due to 
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FIG. 11: (Color online) (a) Contour plot of the station- 
ary charge density of the current-carrying states, in units 
of 10^^ cm~"^. The upper (lower) panel shows the single- 
particle electron density in the QWRs when the electron in 
the quantum dot occupies the bonding (anti-bonding) state. 
The quantum dot tunneling coupling amounts to t = 10 fieV. 
(b) Same type of figure, but with the Coulomb interaction be- 
tween the electron in the wires and the electron in the quan- 
tum dot switched off. In case (b), the two possible paths 
of the scattered electron can interfere with one another. In 
case (a), however, the electron propagating in the lower wire 
changes the double-dot state and looses its capability to in- 
terfere with the other path, thus reducing the visibility (see 
text). 



their charge in this case, the path of the interferometer 
particle is not measured by the double-dot. 

Figures [TT] (a) and (b) show the electron density 
Eq. (P7)) in the entangled Mach-Zehnder double quan- 
tum dot device. To be precise, the figures depict the 
projection of the probability densities of the current- 
carrying two-particle states onto the eigenstates of the 
isolated double quantum dot for a tunneling coupling of 
t — 10 /leV and vanishing phase gate voltage V^sg(P)- 
In particular. Fig. [TT] (a) shows the interacting case, 
whereas Fig. [TT] (b) shows the non- interacting case where 
the electron-electron interaction between the electrons in 
the Mach-Zehnder interferometer and the quantum dot 
electron has been switched off. Thus, the two figures 
illustrate the effect of the Coulomb interaction on the 
charge density. In each of the two figures, the upper 
panel corresponds to the ground state, the lower panel 
corresponds to the excited state of the double quantum 
dot, respectively. In the interacting case (a), we find a 
finite probability of the electron in the double quantum 
dot to be in the excited state. This has to be compared 
to the case of vanishing interaction (b) as well as to the 
case of zero tunneling coupling. In each of the latter 
cases, the probability of the electron in the double quan- 
tum dot to be in the excited state is zero. This is shown 
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in Fig. [Tl](b). Note that for cases (a) and (b), the total 
probabiUty densities that are given by the sum of the two 
projected probabihty densities are almost identical near 
the left contacts, but differ significantly near the right 
contacts. In (a), the total probability density near the 
upper right contact is only partially extinct by interfer- 
ence, whereas in (b), the total probability density van- 
ishes near the upper right contact. This indicates that 
the interaction between the Mach-Zehnder interferometer 
and the double quantum dot destroys the self-interference 
of the QWR-electron wave function. 



D. Discussion: Visibility and von Neumann 
entropy of the entangled Mach-Zehnder double 
quantum dot device 

The results that we have obtained for the entangled 
Mach-Zehnder double quantum dot device are indeed 
the result of a two-qubit quantum operation. In order 
to show this, we apply the analytical model developed 
in Sec. IIVI to our two-qubit device. Figure [12] shows 
the transmission probabilities Ti and T2 , computed from 
Eqs. (|47)) and ([48]), for two different values of the quan- 
tum dots tunneling. Indeed, the comparison with Fig. [9] 
shows that the model is able to reproduce the basic 
physics of the device and to yield transmission proba- 
bilities Ti and T2 that qualitatively agree with the re- 
sults obtained by the numerical simulations described in 
Sees, [n] and IIIII In particular. Figs. [5] and [T^ show a 
similar dependence of the visibility on the tunneling cou- 
pling. The reason for the suppression of the visibility is 
the entanglement of the two qubits. In fact, from the 
point of view of the individual qubit, its entanglement 
with the other one is nothing else but decoherence and 
this decoherence partially suppresses the interference of 
the electron wave function in the two QWRs. In other 
words, the reduced density matrix of each of the qubits, 
obtained by tracing out the degrees of freedom of the 
other qubit, represents a mixed state if and only if the 
two qubits are entangled. Thus, the observed suppression 
of the visibility is both a direct measure of the degree of 
entanglement between two qubits and of the decoherence 
undergone by the QWRs qubit. 

We stress that our approach does not include other 
sources of decoherence as, for example, electron-phonon 
interaction or charge fluctuations in the metallic split 
gates that define the structure. Therefore, our reduction 
of the visibility is obviously ascribed to the quantum en- 
tanglement of the two qubits. On the other hand, in real 
experiments the visibility will be reduced by the coupling 
of the interferometer with any external degree of freedom 
and not only with the other qubit. What makes our re- 
sults functional to experiments is the estimation of the 
visibility decrease against tunable parameters of the re- 
alistic device, as, for example, the tunneling coupling of 
the second qubit. In fact, it is unlikely that the effect 
of the environment will be tuned by the same parame- 




1 2 3 4 5 6 
Phase angle ^ [rad] 

FIG. 12: Transmission functions Ti and T2 as a function of the 
phase angle in rad for the analytical model of the entangled 
Mach-Zehnder double quantum dot device. The full line is 
for a finite tunneling coupling t = 50 in units of /2mL^ , 
whereas the dashed curve is for vanishing tunneling coupling 
that corresponds to no entanglement. 



ters that affect the amount of entanglement between the 
two qubits. The phonon effects on the interferometer, 
for example, are determined solely by the temperature, 
the structure composition and the energy of the travel- 
ing electrons. Thus, our proposed device demonstrates - 
from a modeling perspective - a controlled generation of 
two-qubit entanglement. The predicted behavior of the 
visibility against various system parameters may serve as 
an entanglement witness in experimental realizations of 
the system. 

In the remainder of this section, we show that the re- 
lation between the visibility of our device, as defined in 
Eq. (|50p and depicted in Figs. [TUl and [T^ and the entan- 
glement can be quantified in terms of the von Neumann 
entropy of the reduced density matrix^'"^- that is a mea- 
sure of the degree of entanglement. 

To simplify the derivation, we generally assume that 
back-scattering is negligible in the directional couplers 
and phase gates. This reduces all matrices to size 4x4 
and allows the analytical evaluation of the transmissions 
Ti and T2 in Eqs. ([47]) and ([48]). The density matrix of 
the final state is then given by 

P= H E E^^^yE;iy,\a)'^\Yf{a'r{Y'\''. 

a,a'=0,l Y,Y'=A,B 

(51) 

The amount of entanglement of the final two-qubit state 
is given by the von Neumann entropy of the reduced den- 
sity matrix, namely 



-Trpiod logprcd . 



(52) 



where the reduced density matrix picd is obtained from 
p by tracing out the electron in the double quantum dot. 
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This yields the 2x2 matrix 



Pred 




(53) 



where the so-called coherence Z is given by 

With the eigenvalues of the reduced density matrix, 



(54) 



z\ 



1 



the von Neumann entropy reads 

S = - ^ Pi \ogPi 



(55) 



(56) 



We now show that the eigenvalues p± of piod can be ex- 
pressed in terms of the visibility v, Eq. (|50p. of the inter- 
ference pattern according to 



1 



P± 



(57) 



First, we determine the explicit expressions for the trans- 
mission functions as a function of the angle 4> of the phase 
gate. To this end, we first apply the n/2 rotation gate 
R and the two-qubit gate V onto the initial state I'&in)- 
The probabilities |Ao-,yp of the four two-particle basis 
states |cr)^ |y)° {a = 0,1, Y = A,B) of the resulting 
state are given by 

\Aa,B 
\Aa,A 



icos^ (i(/3±7)), 



(58) 



where -I- and — correspond to cr = and cr = 1, respec- 
tively, and (3 and 7 are angular constants that depend 
on the two-particle interaction V . They describe the lin- 
ear combination of the dot-electron states \B)^ and \A)^ 
within the two-particle wave function. The explicit ex- 
pressions of [3 and 7 do not enter the final result of this 
section. The expressions Eq. ([55]) reflect the fact that the 
probability for the electron in the QWR to be either in 
state |0)^ or in state |1)^ are equal. The next step con- 
sists in applying the phase gate P and the second 7r/2 
rotation gate R. This procedure yields the final state 
probability amplitudes 



Ely = (exp {i(t>) Ao,y + Aiy) , 
El y = (~ ^'^P (*'^) ^o,y + ^i,y) 



(59) 



Note that the gates R and P, being single-qubit trans- 
formations, do not change the degree of entanglement of 



the two electrons. For the transmission probabilities Ti 
and T2 we finally obtain 

Tia = \± \AoM \\Ai,b\ cos (0 + 5b) ^g^^ 
± \Aq^a\ |^i,a|cos {4> + 5a) , 

where the signs ± correspond to Ti and T2, respec- 
tively. The phase angles are defined by 5a = arg(Ao.A) — 
arg(yli.^) and analogously for 5b- By inserting the ex- 
pressions for A„^Y into Eq. (|60p . the transmission prob- 
abilities read 



Ti, - 



1 ± Y cos-^ 13 sin 5 + cos^ 7 cos^ 5 cos (0 -I- (/)o) 

(6i) 

where 5 — ^{5b — 5a) and 0o is a constant phase shift. 
This result shows that the transmission probabilities Ti 
and T2 yield a sinusoidal oscillation as a function of the 
phase angle (j). By comparing the definition Eq. ([50)) of 
the visibility with this result for the transmission prob- 
abilities, we see that the visibility can be written in the 
form 



V — y cos^ (3 sin 5 + cos^ 7 cos^ 5. 



(62) 



T1T2 within 



Now, we evaluate the expression \Z\ 
Eq. ([55|l . and obtain 

\Z\^ - T1T2 = i (cos^ 13 sin^ 5 cos^ 7 cos^ ,5-1) 



1 ( 2 

4^" 



1 



(63) 



Inserting this result into Eq. (j55p . we finally obtain the 
assertion Eq. ([57| . 

This result proves that the von Neumann entropy of 
the reduced density matrix depends monotonously on the 
visibility. A visibility of = 1 corresponds to vanish- 
ing von Neumann entropy (zero degree of entanglement), 
while for a visibility of w = 0, we find a von Neumann 
entropy of S' = 1 (maximum degree of entanglement). 
While the relation Eq. ([57]) between visibility v and von 
Neumann entropy S may not hold rigorously for the re- 
alistic entangled Mach-Zehnder double quantum dot de- 
vice, it is nevertheless interesting to calculate 5* for this 
device as a function of the tunneling coupling, invoking 
the relation of Eq. ([57)l . The result is presented in Fig. [13] 
and shows that the degree of entanglement of the realistic 
device monotonously increases with increasing tunneling 
coupling. 



VI. CONCLUSION 

We have theoretically analyzed semiconductor single- 
and two-qubit quantum gates based on electrostatically 
defined QWRs and quantum dots. We predict the de- 
tailed three-dimensional geometry, material composition. 
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5 10 15 20 
Tunneling coupling [[ieV\ 

FIG. 13: Calculated von Neumann entropy, as defined in the 
text, of the entangled Mach-Zehnder double quantum dot de- 
vice of Sec. |V] C, as a function of the tunneling coupling in 
AteV. 



doping profile, and bias voltage of an all-electric Mach- 
Zehnder interferometer. Our calculations of the elec- 
tronic structure and ballistic transport properties of this 
device show that the proposed device is a fully control- 
lable single-qubit gate for electrons that propagate in 
QWRs. The fabrication of the Mach-Zehnder interfer- 



ometer is within the reach of present-day technology but 
full control of the electronic beam-splitters requires a very 
high fabrication precision in the range of a few nm. 

Based on this all-electric Mach-Zehnder interferome- 
ter, we predict a two-qubit quantum transport device 
where the electrons in the interferometer couple to a 
single-electron double quantum dot by Coulomb inter- 
action. We have calculated the ballistic transport prop- 
erties of the three-dimensional two-qubit device with an 
interacting two-particle quantum transport method and 
designed the device geometry for optimal entanglement. 
By means of an analytical model of the device, we have 
illustrated the qualitative physics of the two-qubit device 
and showed that the visibility is a faithful measure of the 
entanglement. In particular, we have found that the vis- 
ibility can be controlled externally by tuning the tunnel- 
ing coupling between the two quantum dots of the second 
qubit. The device realizes a non-trivial two-qubit gate 
that allows the controlled generation and straightforward 
detection of entanglement from DC current- voltage char- 
acteristics. 
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